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Abstract 

A  low-Reynolds  number  rectilinear  analog  of  the  retreating- 
blade  problem  is  considered  by  computationally  and  exper¬ 
imentally  studying  a  NACA0012  blade  in  span  wise  oscilla¬ 
tion  in  a  free  stream.  Three-dimensional  hybrid  RANS-LES 
simulations  with  spanwise  periodic  boundary  conditions  and 
experimental  flow  visualization  support  the  description  of  ex¬ 
perimental  direct  force  measurements  for  a  wide  range  re¬ 
duced  frequencies  and  advance  ratios,  including  fully  reversed 
flow  conditions.  A  fixed  incidence  of  6  degrees  is  taken  as 
a  nominally  attached-flow  case,  and  agrees  reasonably  well 
with  Isaacs’  theory.  A  fixed  incidence  of  20  degrees  is  taken 
as  a  fully- separated  case,  and  departs  markedly  from  invis- 
cid  theory,  and  even  more  so  from  quasi- steady  approxima¬ 
tion.  Experimental-computational  comparison  shows  a  com¬ 
putational  overprediction  of  lift  relative  to  experimental  re¬ 
sults,  at  moderate  advance  ratios.  Agreement  in  fully  reversed 
flow  is,  however,  quite  good. 


Introduction 

A  major  innovative  thrust  in  rotorcraft  is  the  extension  and 
development  of  both  high  speed  (above  250  knots)  and 
environmentally-friendly  vehicles.  For  conventional  heli¬ 
copters,  the  maximum  speed  has  been  limited  to  advance  ra¬ 
tios  ( ju  =  Uoo/£lR)  at  and  below  0.5,  due  to  the  negative  ef¬ 
fects  associated  with  the  retreating  blade.  Designs  that  incor¬ 
porate  traditional,  compound,  or  coaxial  rotors  (e.g.,  Carter- 
copter  (Ref.[l|  and  the  Sikorsky  X2  Technology  (Ref.  [2]))  may 
apply  slowed  rotational  speeds  to  help  achieve  advance  ratios 
that  approach  and  exceed  1.0. 

The  consequence  of  the  higher  advance  ratios,  exacerbated 
with  slowed  rotation  rates,  is  a  rapid  expansion  of  the  region 
where  reverse  flow  predominates.  In  the  reverse  flow  region, 
the  rotor  blades  operate  with  the  trailing  edge  in  the  relative 
wind  so  that  the  blade  operates  more  as  a  flat  plate  or  bluff 
body  rather  than  an  airfoil.  Reverse  flow  operation  requires 
investigation  and  quantification  to  accurately  capture  these 
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physics  as  they  influence  the  performance  predictions.  Fig- 
ure[l]illustrates  this  increased  region  by  comparing  the  reverse 
flow  region  for  advance  ratios  of  0.5  and  1.0.  The  increase  in 
the  reverse  flow  area  includes  both  inner  and  outer  portions  of 
the  blade,  spanning  a  large  range  of  Reynolds  numbers,  and 
encompassing  both  incompressible  and  compressible  flows. 


Previously,  at  lower  advance  ratios,  these  very  small  re¬ 
verse  flow  regions  were  less  important  in  the  overall  under¬ 
standing  and  prediction  of  the  rotor  behavior.  Johnson  et 
al.  (Ref.  [3]),  Harris  (Ref.  [4])  and  Ormiston  (Ref.  [5])  invigorated 
research  in  reverse  flow  by  highlighting  their  importance  in 
future  vertical  lift  design.  High- advance-ratio  research  on  ro¬ 
tating  systems  has  demonstrated  the  complexity  of  the  flow 
field.  Since  2008,  there  have  been  three  sets  of  high-aspect- 
ratio  wind  tunnel  tests  on  a  three-bladed  NACA0012  rotor 
(Refs.  [6|[7]),  four-bladed  NACA0012  rotor  (Refs.  [8j[9])  and  the 
UH-60  four-bladed  rotor  (Refs.  [Mm  which  have  been  in¬ 
cluded  or  been  accompanied  by  computational  analyses,  for 


example,  but  not  limited  to,  Refs.  9,T2  13 


Beyond  the  rotor  performance  aspects,  at  these  higher  ad¬ 
vance  ratios,  the  reverse  flow  region  can  include  large  regions 
of  radial  flow,  dynamic  stall,  and  surge.  Researchers  have 
appropriately  taken  a  building  block  approach  to  understand 
these  various  mechanisms  in  static  and  non-rotating  dynamic 
situations  prior  to  adding  rotation.  This  approach  has  been 
proven  to  provide  significant  advances  in  understanding  for¬ 
ward  flight  complex  phenomena  such  as  dynamic  stall  (see 
for  example  the  collaborative  efforts  of  the  US,  French  and 
Germans  in  this  area).  If  the  physical  mechanisms  of  each 
are  not  well-understood  for  simpler  flight  conditions  in  re¬ 
verse  flow,  then  incorrect  modeling  or  physical  conclusions 
may  follow.  With  respect  to  radial  or  cross  flow  in  the  rever¬ 


sal  region,  Hodara  and  Smith  (Ref.  14)  and  Bowen-Davies 
and  Chopra  (Ref.  [13])  have  made  advances  to  the  state  of 
the  art  in  modeling  and/or  understanding  the  physics.  Ho¬ 
dara  and  Smith,  for  example,  have  isolated  and  quantified  the 
influence  of  cross  flow  on  airfoil  integrated  quantities  (lift, 
drag,  moment)  in  reverse  flow  and  developed  new  algorithms 
for  comprehensive  codes,  reducing  errors  from  30  %-50  %  to 
less  than  10%.  In  a  recent  collaborative  effort,  the  Univer¬ 
sity  of  Maryland  (UMD)  and  the  Georgia  Institute  of  Tech¬ 
nology  (GIT)  have  partnered  to  quantify,  both  experimen- 
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tally  and  computationally,  static  and  dynamic  stall  in  reverse 
flow  (Ref.  15).  They  have  independently  (Refs.  16  T7])  ob¬ 
tained  similar  static  solutions  and  have  demonstrated  that  a 
reverse  flow  airfoil  can  be  categorized  into  three  distinct  be¬ 
haviors  as  angle  of  attack  increases.  Dynamically  (with  stall), 
the  airfoil  has  five  stages  in  the  reverse  flow  cycle.  Both  static 
and  dynamic  behaviors  in  the  reverse  flow  region  are  different 
than  their  forward  flight  counterparts,  as  described  in  Hodara 
etal.  (Ref.[T5J. 

This  effort  tackles  the  final  major  flow  phenomenon  of  in¬ 
terest  in  reverse  flow,  referred  to  as  surge,  or  the  operation  of 
rotor  blades  in  oscillating  flow  fields.  Contrary  to  dynamic 
stall,  here  the  wing  or  rotor  angle  of  attack  may  remain  either 
constant  or  change,  while  the  free-stream  velocity  oscillates. 
Greenberg  (Ref. [18])  and  Isaacs  (Ref. [19])  pioneered  the  first  ef¬ 
forts  in  this  field  for  forward  flight.  For  forward  flight  config¬ 
urations,  this  topic  has  been  investigated  and  reported  in  mul¬ 
tiple  sources,  the  most  widely  being  Leishman  (Ref.  |20])  and 


Harris  (Ref.  21 ),  which  include  both  theory  and  computation. 
Harris  notes  that  the  addition  of  compressibility  and  viscosity 
(via  high-fidelity  computations)  does  not  significantly  change 
the  magnitude  of  the  lift,  but  does  result  in  a  translation  of  the 
lift  mean  value.  Viscosity  introduces  a  phase  lag  compared  to 
incompressible  theory.  The  physics  associated  with  oscillat¬ 
ing  free-stream  velocity  has  not  been  widely  investigated  for 
reverse  flow  conditions,  although  it  may  be  important  for  new, 
innovative  designs  being  proposed  across  the  engineering  dis¬ 
ciplines,  including  rotorcraft.  Initial  assessment  of  water  tun¬ 
nel  experiments,  theory  (Isaacs/Greenberg),  and  high-fidelity 
computations  are  examined  to  determine  fundamental  differ¬ 
ences  and  areas  where  further  research  are  needed,  building 
on  recent  comparisons  of  theory  and  experiments  (Ref. ! 


Methodology 


Experimental  Setup 

The  U.S.  Air  Force  Research  Laboratory’s  (AFRL)  Horizontal 
Free- surface  Water  Tunnel,  shown  in  Figure [2]  is  fitted  with  a 
three  degree-of-freedom  electric  rig  (Ref.  [24[)~  here  used  only 
in  fore-aft  translation  of  the  test  article  via  the  surge  motor. 
Motion  is  controlled  via  a  Galil  DMC-4040  motion  controller 
from  pre-programmed  scripts,  resulting  in  less  than  0.1mm 
linear  and  0.2°  angular  position  error.  The  test  article,  shown 
in  Figure  [3]  is  a  3D-printed  ABS  plastic  NACA0012  airfoil 
with  a  chord  of  100  mm  and  span  of  457  mm,  leaving  a  nomi¬ 
nally  0.5  mm  wide  gap  between  each  wingtip  and  the  test  sec¬ 
tion  wall.  The  model  is  reinforced  along  the  span  internally 
with  two  6.5  mm  diameter  carbon  fiber  tubes.  The  blade  is 
installed  upside  down  in  the  water  tunnel  so  that  positive  lift 
acts  downwards  in  Figure  [3] 

Flowfield  visualization  was  performed  using  dye- 
illuminated  planar  laser  fluorescence.  A  high  concentration 
of  Rhodamine  6G  in  water  was  injected  at  the  leading  and 
trailing  edges  at  3/4-span.  Dye  flow  rate  was  controlled  by 
a  positive-displacement  pump,  connected  to  a  set  of  0.5  mm 
internal-diameter  rigid  lines  glued  to  the  surface  of  the  plate 


fjV« 


Fig.  1.  Iso-sweep  angles  [deg]  over  the  rotor  disk  at  /i  =  0.5 
and  1.0  (counter-clockwise  rotation).  The  red  line  indi¬ 
cates  the  reverse  flow  region  (Ref.  [23). 

as  documented  by  (Ref.  [24]).  The  dye  was  illuminated  by 
an  Nd:YLF  527  nm  pulsed  laser  sheet  of  2  mm  thickness 
at  50  Hz,  and  images  were  recorded  with  a  PCO  DiMax 
high-speed  camera  through  a  Nikon  PC-E  45  mm  Micro  lens. 
A  Tiffen  orange  21  filter  was  used  to  remove  the  incident  and 
reflected  laser  light,  leaving  only  dye  fluorescence. 

Force  measurements  were  performed  using  a  submergi¬ 
ble  ATI  Nano25  IP68  6-component  strain  gauge  force/torque 
sensor,  sampled  at  100  Hz  with  a  hardware  low-pass  filter  at 
34  Hz,  and  then  it  was  filtered  at  5  times  the  motion  frequency 
using  a  4th-order  Chebyshev  II  low-pass  filter  in  Matlab.  All 
motions  were  repeated  for  one  hundred  cycles;  the  first  three 
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Fig.  2.  AFRL  free-surface  water  tunnel  test  section  with 
planar  motion  mechanism. 


Fig.  3.  NACA0012  airfoil  installed  in  the  AFRL  free- 
surface  water  tunnel. 

cycles  were  not  included  in  the  ensemble  averages  given  here 
to  remove  transients.  In  all  of  the  plots  below,  time  has  been 
nondimensionalized  by  the  motion  period  T,  while  Q  has 
been  nondimensionalized  by  the  free  stream  in  the  water  tun¬ 
nel,  Uoo. 

Blade  Kinematics 

In  the  current  work,  a  rotating  rotor  blade  has  been  modeled 
as  a  two-dimensional  blade  element  in  harmonic  motion.  Ne¬ 
glecting  inflow  and  flapping,  the  relative  flow  over  a  blade  el¬ 
ement  of  a  rotor  in  hover  is  given  by  U  =  Tlr,  where  £1  is  the 
angular  velocity  of  the  rotor  blade  and  r  is  the  radial  location 
of  the  blade  element.  For  a  rotor  in  forward  flight,  the  addition 
of  a  relative  free  stream,  Uoo,  introduces  a  harmonic  oscillation 
and  the  incident  flow  over  a  blade  element  becomes 

U(y)  =  £lr  +  Uoo  sin  \\f  =  CLr-\-  jiTlR sin  y/  (1) 

where  /I  =  Uoo/ £IR  is  the  advance  ratio  and  R  is  the  rotor  ra¬ 
dius.  Equation  [T]  is  plotted  in  Figure  [4j  where  the  azimuthal 
position  of  the  blade  (i.e.,  0°  <  i jf  <  360°)  is  mapped  to  a 
nondimensional  time  scale  t/T.  The  maximum  relative  flow 
Umax  =  Uoo  +  £lr  =  £l(r  +  fdR)  occurs  at  a  blade  azimuthal  an¬ 
gle  of  1 ff  =  90°,  corresponding  to  t/T  =  0.25.  The  minimum 


Fig.  4.  Normalized  velocity  during  one  rotor  cycle. 

relative  flow  Um\n  =  Uoo  —  Tlr  =  Q(r  —  fiR)  occurs  at  a  blade 
azimuthal  angle  of  if/  =  270°,  corresponding  to  t/T  =  0.75. 
The  amplitude  of  the  oscillation  of  the  relative  flow  experi¬ 
enced  by  a  blade  element  is  the  free  stream,  Uoo,  and  the  mean 
velocity  at  radial  potion  r  is  U  =  £lr. 


To  model  the  rotational  motion  as  rectilinear  motion,  the 
amplitude  of  the  velocity  oscillation  is  nondimensionalized  by 
the  local  mean  velocity  to  define  a  blade-element  advance  ra¬ 
tio 


JL 

r/R 


(2) 


where  R  is  the  rotor  radius,  /I  =  U^/TIR  is  the  advance  ratio 
of  the  rotor,  and  A  is  a  parameter  that  characterizes  the  local 
(radial)  magnitude  of  reverse  flow.  When  r  =  R,  A  =  /I,  and 
when  r  <  R,  A  >  ju.  Blade  elements  inboard  of  the  blade  tip 
experience  more  extreme  reverse  flow  than  does  the  blade  tip. 
For  example,  if  jl  =  1,  then  A  =  /i  =  1  at  the  blade  tip  (r/R  = 
1).  At  the  mid- span,  however,  r/R  =  0.5  and  A  =  2/i.  In  this 
case,  the  most  severe  reverse  flow  experienced  at  the  tip  of  the 
rotor  blade  is  a  relative  velocity  of  zero,  while  the  mid-span 
experiences  a  reverse  flow  that  is  twice  the  free  stream.  A  also 
characterizes  the  radial  extent  of  reverse  flow  for  a  given  value 
of  /i,  as  illustrated  in  Fig.  [5]  For  example,  if  jd  =  0.5,  A  >  1 
for  r/R  <  0.5.  Regions  of  the  blade  inboard  of  r/R  =  0.5  will 
experience  reverse  flow. 


To  produce  a  harmonic  oscillation  in  rectilinear  motion 
analogous  to  that  experienced  by  the  rotating  blade  element, 
the  blade  mounted  in  the  water  tunnel  was  driven  forward  and 
aft  in  a  steady  free  stream  of  speed  Uoo  with  a  velocity  profile 
of 


u(t/T)  =  Uoo  A  sin  (cot/T)  (3) 


where  (O  is  the  angular  frequency  of  the  nondimensional  time 
scale,  t/T.  Note  that  u'(t)  is  the  speed  at  which  the  model  is 
driven  and  varies  from  —Uoo  to  -\-Uoo  for  A  =  1.  In  the  wa¬ 
ter  tunnel,  the  harmonic  oscillation  driven  by  the  surge  mo¬ 
tor  represents  the  forward-flight  contribution  to  velocity.  The 
unsteady  flow  incident  to  the  blade  element  is  computed  by 
adding  the  component  from  the  rotation  of  the  rotor  such  that 


u(t/T )  =  0)r  +  f/ooA  sm((Qt/T)  =  £/«>[  1  + A  sin  (cot/T)]  (4) 
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Table  1.  Test  matrix.  Force  histories  are  available  for  all 
of  the  cases  listed.  Dye  flow  visualization  and  numerical 
simulations  at  0  =  20°  are  available  for  the  cases  in  italics. 


Reduced  frequency 

k  =  nfc/Uoo 

Blade-element  advance  ratio 

A  =  (2  nAf/Uco)  =  2 Ak/c 

0.0952 

0.57, 0.71,0.86,  1.00 

0.1111 

0.67,0.83,  1.00,  1.17 

0.1333 

0.80,  1.00,  1.20,  1.40 

0.1667 

1.00,  1.25,  1.50,  1.75 

0.2222 

1.33,  1.67,  2.00,  1.75 

Note  that  in  the  water  tunnel,  the  free  stream  Uoo  represents 
the  velocity  due  to  blade  rotation,  and  so  Uoo  =  cor.  If  A  =  1 
as  in  the  previous  example,  the  blade  element  experiences  a 
relative  flow  ranging  from  0  to  2Uoo. 

The  position  of  the  blade  element  in  the  driven  rectilinear 
harmonic  motion,  s(t)  is  given  by  the  integral  of  the  velocity 
profile  in  Eq.  [3j 

s(t/T )  =  ~ ~“~cos  (cot/T)  (5) 


The  physical  frequency  of  the  blade  position,  /  =  co/(2n), 
can  be  written  in  terms  of  the  reduced  frequency,  k  =  (Oc/2Uoo , 
as 


/  =  kUoo  /  7l  c 

and  the  amplitude  of  the  blade  motion  as 

Uook  _  Ac 

co  2k 


A  = 


(6) 

(7) 


Thus,  for  a  fixed  blade  chord  c  and  free  stream  Uoo,  the  param¬ 
eter  space  is  defined  by  the  blade-element  advance  ratio  A  and 
the  reduced  frequency  k. 

Experiments  were  performed  in  water  at  a  constant 
Reynolds  number  of  Re  =  Uooc/v  =  40000,  thereby  fixing  the 
value  of  Uoo.  The  amplitude  of  the  blade  motion  was  limited 
by  the  0.6  m  stroke  length  of  the  surge  motor.  The  reduced 
frequencies  tested,  0.095  <k<  0.222,  were  limited  by  the  ca¬ 
pabilities  of  the  force  balance.  The  experimental  parameter 
space  in  k  and  A  is  given  in  Table  [I]  and  illustrated  in  Fig.  [5] 
Experiments  were  run  with  the  NACA0012  blade  mounted  at 
an  angle  of  incidence  of  6  =  6°  and  20°  to  assess  behavior  in 
nominally  attached  and  separately  conditions,  respectively. 


Numerical  Methods 

A  numerical  solver  has  been  developed  at  GIT  to  study  un¬ 
steady  problems  such  as  dynamic  stall  and  rotor  blades  in 
reverse  flow.  This  new  platform  was  specifically  designed 
to  rapidly  implement  and  evaluate  state-of-the-art  algorithms, 
including  turbulence  and  transition  approaches.  The  three- 
dimensional  compressible  governing  equations  are  solved  in 
a  time-accurate  manner  using  a  cell-centered  finite  volume 
approach  on  structured  grid  topologies.  The  spatial  recon¬ 
struction  is  carried  out  using  Van  Leer’s  MUSCL  scheme, 
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Fig.  5.  Text  Matrix  in  the  k  —  A  parameter  space,  with 
symbols  colored  to  represent  C/min  for  a  blade  at  6  =  6°. 
Symbols  represent  points  in  the  parameter  space  at  which 
force  data  was  acquired  for  both  6  =  6°  and  20°.  Sym¬ 
bols  marked  with  an  6x’  represent  cases  for  which  dye 
flow  visualization  and  numerical  simulations  were  also 
performed  (for  6  =  20°). 

leading  to  second-order  accuracy.  The  convective  fluxes  are 
computed  using  Roe’s  flux  difference  splitting  scheme,  while 
the  viscous  fluxes  are  obtained  from  second-order  central  dif¬ 
ferences.  The  solution  is  marched  in  time  using  the  im¬ 
plicit  LU-SSOR  scheme  combined  with  the  method  of  Gear 
to  achieve  second-order  temporal  accuracy.  The  number  of 
Newton  sub-iterations  at  each  physical  time  step  is  chosen  to 
ensure  a  proper  convergence  of  the  residuals.  A  large  num¬ 
ber  of  boundary  conditions  has  been  implemented,  including 
inviscid/viscous  walls,  non-reflecting  inlets/outlets  based  on 
Riemann  invariants,  etc.  The  solver  is  fully  parallel  (MPI)  and 
includes  rigid  body  rotation/translation  using  the  Arbitrary 
Lagrangian  Eulerian  approach  (Ref.  [25]).  In  the  present  study, 
the  turbulent  equations  are  solved  in  a  loosely  coupled  man¬ 
ner  using  the  hybrid  Reynolds -Averaged  Navier  Stokes  -  large 
eddy  simulation  (HRLES)  model  of  Sanchez-Rocha  (Ref.[26|. 
Validation  of  the  code  has  been  verified  over  a  large  number 
of  test  cases,  ranging  from  turbulent  circular  cylinders  to  os¬ 
cillating  wings  in  reverse  flow,  several  of  which  are  illustrated 
in  Hodara  et  al.  (Ref.  m 

In  the  current  work,  numerical  simulations  were  performed 
on  a  semi-infinite  NACA0012  blade  at  an  angle  6  =  20°  with 
respect  to  the  free  stream.  Grid  motion  reproduced  selected 
cases  from  the  experimental  test  matrix,  each  with  an  oscilla¬ 
tion  amplitude  A/c  =  A  /(2k)  =  3.  Maintaining  this  constant 
amplitude,  three  cases  of  different  reduced  frequency  were 
evaluated:  k  =  0.167,  0.133  and  0.095.  Fixing  the  amplitude 
of  the  motion  and  the  reduced  frequency,  A  becomes  1.0,  0.8, 
and  0.57,  respectively,  for  the  three  cases.  Figure [5] illustrates 
where  these  cases  lie  in  the  parameter  space  explored  in  the 
current  work.  The  A  =  1.0  case  is  of  particular  interest  be¬ 
cause  the  relative  velocity  in  a  blade-element- fixed  reference 
frame  reaches  zero  at  t/T  =  0.75,  corresponding  to  the  onset 
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t/T  =  0.25 
u  =  1-8  Un 


(a)  Lift  histories 


t/T  =  0.75 
u  =  0.20  U0 


(b)  Dye  flow  visualization 


Fig.  6.  Experimental  results  for  a  NACA0012  at  0  =  6°.  Dye  flow  visualization  is  shown  for  A0.1333  =  0.8. 


of  reverse  flow. 


Computational  Grid 


The  experimental  test  conditions  discussed  previously 
were  reproduced  in  the  numerical  simulations.  The  Reynolds 
number  based  on  the  airfoil  chord  is  Rec  =  40000  and  the 
free-stream  Mach  number  is  Moo  =  0.1.  This  Mach  number 
does  not  correspond  to  the  experimental  value  (Moo  <0.001) 
due  to  the  density-based  nature  of  the  solver  and  the  current 
absence  of  a  low  Mach  number  preconditioner.  The  Reynolds 
number  and  unsteady  characteristics  were  scaled  to  reproduce 
the  experiment  given  the  change  in  Mach  number.  The  change 
in  density  was  found  to  remain  consistantly  below  4  %  when 
the  dynamic  pressure  was  maximum  (case  Aq.167  =  1-0  at 
t/T  =  0.25).  It  is  therefore  expected  that  compressible  effects 
will  be  negligible  in  the  present  study  based  on  aerodynamic 
theory. 

The  physical  time  step  was  maintained  for  all  three  cases 
evaluated  numerically,  with  a  flow  particle  taking  approxi¬ 
mately  1250  time  steps  to  travel  from  the  leading  edge  to  the 
trailing  edge  of  the  airfoil.  Based  on  the  oscillation  periods  of 
the  three  cases,  the  number  of  time  steps  per  cycle  is  25  000, 
30000  and  40000  for  A  =  1.0,  0.8  and  0.57  respectively.  All 
simulations  included  the  HRLES  turbulence  closure,  and  the 
results  were  phase-averaged  over  12,  7,  and  5  cycles  for  the 
three  respective  cases.  The  number  of  cycles  becomes  small 
for  lower  values  of  A,  but  smaller  A  coincides  with  smaller 
fluctuations  in  forces  and  moments  from  cycle  to  cycle,  which 
mitigates  the  need  for  ensemble  averaging  over  large  numbers 
of  cycles.  Due  to  the  fine  time  steps  employed  in  the  simula¬ 
tions,  5  Newton  subiterations  were  sufficient  for  reducing  the 
L°°  norm  of  the  residuals  by  two  orders  of  magnitude  between 
each  physical  time  step.  A  temporal  analysis  concluded  that 
increasing  the  number  of  subiterations  further  had  a  negligible 
impact  on  the  results,  similar  to  prior  conclusions  by  Liggett 
and  Smith  (Ref.  [27]). 


A  rigorous  grid  refinement  study  was  conducted  to  identify 
the  optimal  mesh  for  the  present  configuration;  a  family  of 
O-grids  was  generated  following  the  guidelines  of  Spalart 
(Ref.  |28)  and  Smith  et  al.  (Ref.  [16]).  The  spacing  at  the  wall 
was  chosen  to  ensure  that  y+  <  1  in  the  first  cell,  with  approxi¬ 
mately  50  points  resolving  the  boundary  layer.  The  number  of 
points  in  the  circumferential  direction  was  increased  from  499 
for  the  coarse  grid  to  999  for  the  fine  grid,  while  the  number 
of  grid  points  in  the  wall  normal  direction  was  increased  from 
161  to  201.  The  computational  domain  extends  approximately 
50  chord-lengths  from  the  blade  section  in  all  directions.  The 
radial  or  spanwise  dimension  of  the  domain  must  be  carefully 
chosen  to  ensure  that  the  periodic  boundary  conditions  do  not 
impact  the  flow  field  in  three-dimensional  simulations.  Most 
researchers  (Refs.|29}]3T])  have  used  a  width  of  1  x  c  to  simu¬ 
late  the  flow  past  a  semi-infinite  wing  in  post-stall  conditions, 
while  Smith  (Ref.  |T6|  used  a  larger  width  of  2  x  c.  In  the 
present  study,  the  coarsest  grid  has  a  width  of  1  x  c,  which 
is  extended  to  2  x  c  for  the  finest  mesh.  In  all  cases,  the  ra¬ 
dial  aspect  ratio  is  held  constant,  with  31  points  per  chord 
length  (z+  cx  200  at  the  Reynolds  number  considered).  This 
gridding  approach  was  found  to  generate  accurate  (when  cor¬ 
related  with  water  and  wind  tunnel  experimental  data)  static 
and  dynamic  simulations  in  forward  and  reverse  flow  condi¬ 
tions,  as  demonstrated  by  Smith  et  al.  (Ref.[T6|)  and  Hodara  et 
al.  (Ref.|T3). 


Several  meshes  were  evaluated  in  the  simulation  of  the 
flow  past  a  semi-infinite  static  NACA0012  wing  at  Rec  = 
1.1  x  105,  deep  into  the  reverse  stall  regime  (a  =  150°).  This 
challenging  configuration  was  computed  using  the  HRLES 
turbulence  closure,  with  approximately  5000  time  steps  re¬ 
solving  each  vortex  shedding  cycle.  Lift  and  drag  coefficients, 
as  well  as  Strouhal  number  were  captured  within  3  %  of  one 
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another  for  the  two  meshes  identified  here.  The  solution  ob¬ 
tained  with  the  coarser  mesh  (499  x  31  x  161  over  a  single 
chord-width  span  dimension)  is  considered  sufficient  and  ap¬ 
plied  in  the  present  study. 

Results  and  Discussion 

Experiments  were  performed  on  a  NACA0012  blade  mounted 
at  6  =  6°  and  20°  incidence  over  a  wide  k  —  A  parameter 
space,  given  in  Table  [T]  and  Figure  [5]  Numerical  simulations 
were  performed  on  a  semi-infinite  NACA0012  wing  at  an  in¬ 
cidence  angle  of  0  =  20°.  Of  particular  interest  is  the  transi¬ 
tion  from  forward  to  reverse  flow  that  exists  at  A  =  1.0,  where 
the  relative  velocity  in  a  blade-element-fixed  reference  frame 
reaches  zero  at  t/T  =  0.75.  The  two  incidence  angles  were 
selected  to  provide  one  case  expected  to  be  free  of  significant 
flow  separation  suitable  for  comparison  with  inviscid  theory, 
and  one  case  more  representative  of  the  separated  flow  likely 
to  exist  in  the  reverse  flow  region  of  a  rotor  disk. 


Results  at  Low  Incidence 


Figure  |6(a)  gives  the  experimentally  measured  phase- 
averaged  lift  coefficient  as  a  function  of  nondimensional  time 
for  the  NACA0012  blade  at  0  =  6°.  The  lift  coefficient  is 
based  on  the  free  stream  condition  independent  of  the  oscil¬ 
latory  motion,  i.  e.,  Cl  =  2L/ (pU^S)  where  Uq  is  the  tunnel 
freestream  and  S  is  the  wing  area.  The  lift  coefficient  history 
is  shown  in  the  figure  as  the  instantaneous  unsteady  lift  coeffi¬ 
cient  Cl,  normalized  by  the  static  lift  coefficient  as  measured 
at  G  =  6°,  Cl0  =0.57.  Because  both  the  unsteady  and  static 
lift  coefficients  are  based  on  the  free  stream  velocity  of  the  wa¬ 
ter  tunnel,  Cl/Cl0  =L/Lq.  Note  that  in  the  following  text  and 
figures,  the  subscripts  on  A  are  the  value  of  k  for  that  case. 
For  example,  A0.222  =  1-33  refers  to  a  case  where  A  =  1.33 
and  k  =  0.222. 


Force  curves  are  shown  in  Figure  |6(a)| for  each  of  the  cases 
evaluated,  with  the  one  corresponding  to  the  dye  flow  visu¬ 
alization  in  Figure  |6(b)|  highlighted  in  green.  In  each  case, 
lift  rises  as  the  blade  surges  into  the  free  stream  and  the  rela¬ 
tive  flow  increases.  The  maximum  relative  velocity  occurs  at 
t/T  =  0.25,  and  the  maximum  lift  occurs  just  past  this  point, 
at  t/T  =  0.28.  Lift  decreases  as  the  blade  retreats  and,  for 
most  cases,  remains  near  zero  for  a  substantial  portion  of  the 
retreating  phase  of  the  motion  cycle.  Dye  flow  visualization 
of  the  A0.1333  =0.8  case  is  shown  in  Figure [6(b)] at  t/T  =  0.25 
and  0.75.  At  t/T  =  0.25,  the  blade  is  advancing  into  the 
free  stream,  while  the  blade  is  retreating  at  t/T  =  0.75.  As 
A  <  1.0,  the  blade  section  does  not  experience  reverse  flow. 
The  relative  flow  vectors,  u ,  are  shown  to  scale  in  the  figure. 
At  this  incidence  angle,  flow  is  attached  to  the  blade  during 
the  advancing  phase,  though  incipient  separation  is  observed 
at  t/T  =  0.75,  where  the  relative  flow  velocity  is  lowest. 

In  a  few  cases  shown  in  Figure  |6(a)|  lift  becomes  nega¬ 
tive  during  the  retreating  phase  (See  also  Fig.  [5]).  The  largest 
negative  lift  is  produced  by  the  A0.222  =  2.0  case,  followed 


by  A0.222  =  1.67  and  Ao.167  =  1-75.  A  small  amount  of  neg¬ 
ative  lift  is  produced  by  the  A0.222  =  1-33  and  Ao.i67  =  1.5 
cases.  Negative  lift  was  only  observed  for  cases  where  A  is 
significantly  greater  than  1 .  Here,  the  relative  flow  direction 
reverses  and  the  blade  experiences  flow  from  the  geometric 
trailing  edge  to  the  geometric  leading  edge  (i.e.,  right  to  left 
in  Figure  [6(b)]).  In  this  orientation,  a  blade  that  was  at  a  posi¬ 
tive  angle  of  attack  is  now  at  a  negative  angle  of  attack  and  so 
the  direction  of  lift  reverses. 

Results  at  High  Incidence 

Numerical  Flow  Visualization  Numerical  simulations  at 
0  =  20°  provide  instantaneous  flow  fields  throughout  the 
blade  motion,  demonstrating  the  evolution  of  the  flow  struc¬ 
ture.  Two-dimensional  slices  of  the  flow  field  were  ex¬ 
tracted  at  various  times  during  the  cycle  to  investigate  the  flow 
physics  of  the  various  cases.  These  slices  were  taken  dur¬ 
ing  a  single  (last)  cycle  at  the  midspan  location  and  have  not 
been  ensemble-averaged.  Contours  of  vorticity  magnitude  are 
shown  in  Figures  [7]  and  [8]  at  eight  instances  during  the  blade 
motionfor  three  cases.  Note  that  since  A/c  =  A/(2&)  =  3  is 
constant  for  the  three  cases  shown  here,  the  reduced  frequency 
of  the  motion  increases  with  A . 

The  beginning  of  the  cycle,  t/T  =  0,  is  defined  as  the  in¬ 
stant  when  the  blade  section  is  the  furthest  downstream  and 
the  oscillatory  velocity  is  equal  to  zero  (see  Fig.  [4]).  The  rela¬ 
tive  flow  velocity  at  this  instant  is  thus  equal  to  the  free  stream. 
It  can  be  seen  in  row  1  of  Figure  [7]  that  the  flow  structure  at 
this  point  in  time  differs  significantly  from  that  over  a  static 
wing,  and  varies  with  reduced  frequency  of  the  motion.  In 
the  case  where  A  =  0.57,  the  flow  is  separated  from  near  the 
leading  edge,  much  like  a  static  airfoil  at  such  high  incidence. 
As  the  reduced  frequency  increases,  however,  dynamic  effects 
become  more  pronounced  and  the  flow  remains  attached  over 
a  large  portion  of  the  suction  side  of  the  blade.  The  physi¬ 
cal  explanation  for  this  boundary  layer  reattachment  becomes 
apparent  at  the  end  of  the  cycle.  Moving  forward  in  time, 
the  blade  accelerates  into  the  free  stream  and  reaches  its  max¬ 
imum  relative  velocity  when  t/T  =  0.25  (row  3  of  Fig.  [7]). 
At  this  time,  the  flow  features  are  similar  for  all  three  cases, 
though  the  shear  layer  appears  to  be  less  stable  at  A  =  0.8  than 
in  the  other  cases.  In  all  three  cases,  the  boundary  layer  sepa¬ 
rates  close  to  the  leading  edge  and  a  large  region  of  separated 
flow  and  thus  low  pressure  exists  on  the  suction  side  of  the  air¬ 
foil.  These  effects,  combined  with  the  large  dynamic  pressure 
on  the  pressure  side,  yield  lift  forces  up  to  8  times  the  static 
value  (for  A  =  1.0),  as  shown  in  Figure [l2| 

In  the  second  half  of  the  cycle,  the  blade  decelerates  to 
t/T  =  0.5,  where  the  oscillation  velocity  approaches  zero  and 
the  freestream  speed  is  recovered.  The  flow  fields  at  this  point 
are  given  in  row  1  of  Figure  [8]  The  shear  layer  originating 
at  the  leading  edge  extends  away  from  the  blade  surface  at  an 
angle  dependent  on  the  blade-element  advance  ratio  A.  The 
variation  in  the  shear  layer  angle  is  barely  perceptible  in  these 
images  at  t/T  =  0.5,  but  becomes  much  more  significant  as 
the  blade  section  begins  its  retreating  phase,  shown  in  row  2  of 
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t/T  =  0.375  t/T  =  0.25  t/T=  0.125  t/T  =  0.0 


Fig.  7.  Numerically-generated  instantaneous  vorticity  magnitude  contours  on  a  NACA0012  blade  section  at  Re  =  40000, 
0  =  20°  for  0  <  t/T  <  0.5.  The  three  cases  shown  here  are  for  oscillations  defined  by  A0.095  =  0.57,  A0.133  =  0.8,  and 
^0.167  =  1-0. 
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t/T  =  0.875  t/T  =  0.75  t/T  =0.625  t/T  =  0.5 


Fig.  8.  Numerically-generated  instantaneous  vorticity  magnitude  contours  on  a  NACA0012  blade  section  at  Re  =  40000, 
0  =  20°  for  0.5  <  t/T  <  1.0.  The  three  cases  shown  here  are  for  oscillations  defined  by  A0.095  =  0.57,  A0.133  =  0.8,  and 
^0.167  =  1-0. 
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Figure [8]  At  t/T  =  0.625  the  leading  edge  shear  layer  breaks 
down  for  both  A  =0.80  and  1.0,  and  a  low  velocity  turbulent 
region  appears  on  the  suction  side  of  the  airfoil.  This  breakup 
occurs  earlier  for  A  =0.8  than  for  A  =  1.0  due  to  the  more 
unstable  shear  layer  in  the  former  case,  and  may  be  related  an 
effect  of  reduced  frequency.  Remember  that  in  the  current  set 
of  experiments,  the  amplitude  of  the  blade  motion  was  fixed, 
so  both  the  blade-element  advance  ratio  A  and  reduced  fre¬ 
quency  k  vary  together. 

At  t/T  =  0.75  (row  3  in  Fig.  [8]),  the  shear  layer  extending 
downstream  from  the  leading  edge  persists  at  A  =  0.57.  At 
higher  reduced  frequencies,  the  turbulent  wake  front  moves 
slightly  upstream  of  the  airfoil.  This  is  initially  counter¬ 
intuitive  since  the  relative  velocity  in  the  body  frame  does 
reach  zero  at  A  =  1  and  t/T  =  0.75,  but  does  not  become 
negative  at  any  time.  However,  note  that  the  blade  section 
is  retreating  and  so  flow  structures  that  advect  more  slowly 
than  the  freestream  will  appear  to  move  upstream  relative  to 
the  blade  section.  The  numerical  flow  visualization  comports 
well  with  experimental  flow  visualization  (see  Fig.[T0|). 

Finally,  at  t/T  =  0.875  (shown  in  row  4  in  Fig.  [8]),  the 
boundary  layer  begins  to  reattach  at  the  leading  edge  of  the 
blade  as  the  relative  flow  velocity  increases  once  again.  This 
is  most  evident  in  the  A  =  0.57  case  where  the  incident  flow 
has  remained  of  the  largest  magnitude.  In  this  case,  the  shear 
layer  separation  angle  remained  relatively  small  throughout 
the  cycle  and  the  turbulent  wake  did  not  reach  the  leading  edge 
of  the  airfoil.  Throughout  the  cycle,  structure  of  the  flow  field 
for  this  case  remained  more  similar  to  that  of  a  static  airfoil 
than  did  those  at  higher  blade-element  advance  ratios.  At  low 
A,  variations  in  the  incident  flow  velocity  are  small  enough 
that  the  overall  flow  structure  is  not  fundamentally  changed  as 
the  blade  advances  and  retreats  in  the  free  stream.  At  higher 
A,  larger  variations  in  incident  flow  velocity  have  resulted  in 
very  low  flow  velocities  and  near- stagnant  flow  as  the  blade 
retreats.  When  the  blade  then  reverses  direction  and  begins 
the  cycle  again,  there  is  a  small  window  of  time  in  which  an 
attached  boundary  layer  forms  (over  at  least  part  of  the  blade 
section)  in  this  near- stagnant  flow  before  flow  separation  once 
again  occurs. 

The  case  of  A  =  0.8  serves  to  further  illustrate  how  the  flow 
behaves  between  large  separation  (A  =  0.57)  and  attached 
flow  (A  =  1.0)  at  time  t/T  =  0.0.  As  the  cycle  progresses, 
the  instability  that  lies  at  A  =  0.8,  between  the  two  starting 
extremes,  becomes  obvious.  As  the  cycle  approaches  its  mid¬ 
point  (t/T  =  0.5),  the  reverse  flow  from  the  leading  edge  shear 
layer  and  trailing  edge  vortex  interaction  has  penetrated  clos¬ 
est  (of  the  three  A  configurations)  to  the  leading  edge.  For 
A  =  0.8,  the  leading  edge  shear  layer  angle  maintains  a  trail¬ 
ing  edge  vortex  that  is  stronger  and  closer  to  the  airfoil  trail¬ 
ing  edge.  By  t/T  =  0.625,  the  middle  A  =  0.8  case  is  striking 
with  its  unstable  shear  layer.  It  then  breaks  down  differently 
than  at  advance  ratios  greater  and  less  than  this  value.  This 
instability  can  have  significant  impact  on  the  behavior  of  the 
rotor  blade  when  combined  with  angular  motion,  and  could 
partially  explain  why  some  airfoil  and  rotor  configurations  in 


Fig.  9.  Experimental  lift  histories  for  a  NACA0012  blade 
section  at  0  =  20°. 

the  literature  appear  to  be  more  sensitive  during  dynamic  stall 
than  others. 

Experimental  Results  Experimental  results  at  0  =  20°  are 
given  in  Figures  |6(a)|  and |T0|  The  force  histories  have  a  shape 
similar  to  that  observed  at  0  =  6°,  with  lift  rising  as  the  blade 
surges  into  the  freestream  and  dropping  to  near  zero  as  the 
blade  retreats.  At  this  high  angle  of  attack,  however,  maxi¬ 
mum  lift  no  longer  lags  the  maximum  dynamic  pressure.  The 
maximum  lift  recorded  during  the  experiments  appears  to  lead 
the  blade  motion,  and  more  so  as  A  increases.  While  counter¬ 
intuitive,  similar  phase-lead  behavior  has  been  previously  ob- 
served  at  0.  \<k<  0.6  for  a  NACA0009  at  A  =  0. 1  (Ref.  [22). 

Experimental  flow  visualizations  are  available  for  the  cur¬ 
rent  experiments  at  t/T  =  0.25,  given  in  Figure  [To|).  These 
images  exhibit  many  of  the  same  flow  features  observed  in 
the  numerical  simulations.  The  angle  of  flow  separation  evi¬ 
dent  in  the  dye  flow  visualization  at  t/T  =  0.75  is  drastically 
different  for  the  three  cases  as  was  observed  in  the  numer¬ 
ical  results.  It  should  be  noted  that  the  large  vortex  visible 
upstream  of  the  leading  edge  in  the  dye  flow  visualization  at 
A  =  1  was  also  predicted  in  the  numerical  simulations,  albeit 
at  a  later  time  of  t/T  «  0.87.  This  vortex  is  an  artifact  of 
the  leading-edge  shear  layer  rolling  up  as  the  blade  retreats. 
Despite  obvious  differences  in  the  structure  of  the  separated 
flow  across  the  three  cases  shown  here,  the  flow  is  fully  sep¬ 
arated  in  each  case  and  the  blade  produces  near-zero  lift  at 
t/T  =  0.15. 

Figure  El  gives  the  lift  histories  for  four  different  cases 
in  which  A  =  1  is  fixed  and  the  reduced  frequency  k  (and  thus 
the  amplitude  of  the  motion,  A/c)  is  varied.  The  curves  shown 
here  are  a  subset  of  those  given  in  Figure  |6(a)|  Although  the 
maximum  lift  achieved  in  these  cases  varies  in  magnitude,  the 
rate  at  which  lift  builds  up  over  0  <  t/T  <0.1  and  falls  off 
over  0.25  <  t/T  <0.5  is  similar.  In  the  second  half  of  the 
cycle,  0.5  <  t/T  <  1,  the  magnitude  of  lift  agrees  as  well. 
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Fig.  10.  Experimental  dye  flow  visualization  for  a  NACA0012  blade  section  at  0  =  20°. 


Fig.  11.  Lift  histories  for  a  fixed  A  =  1.0  and  varying  k. 


It  should  be  noted  that  the  normalized  apparent  mass  contri¬ 
bution  to  lift,  La/Lo ,  is  proportional  to  the  product  of  blade- 
element  advance  ratio  A  and  reduced  frequency  k.  In  the  cases 
presented  in  this  paper,  apparent  mass  is  a  negligible  contribu¬ 
tor.  The  difference  in  peak  normalized  lift  t/T  =  0.25  in  Fig¬ 
ure  11  for  varying  reduced  frequency  k  at  constant  A  =  1.0, 


must  therefore  come  from  a  difference  in  circulation  and  vor¬ 
tex  advection  speed  from  the  preceding  retreating  phase. 


Comparison  of  Integrated  Forces 


The  lift  histories  as  computed  numerically,  measured  experi¬ 
mentally,  and  predicted  by  inviscid  theory  are  summarized  in 
Figure [12]  In  this  figure,  all  of  the  lift  histories  are  normalized 
as  described  previously,  i.e.,  by  the  viscous  static  lift,  Lq.  For 
all  cases  the  broadest  variations  and  the  most  interesting  phe¬ 
nomena  in  force-histories  occur  on  the  advancing  semi-stroke, 
especially  between  t/T  =  0  and  t/T  ~  0.3. 


Inviscid  theory  does  a  reasonable  job  of  predicting  un¬ 
steady  lift  at  0  =  6°  (Fig.  12(a)|),  though  lift  is  overpredicted 
during  the  first  half  of  the  stroke.  The  theoretical  model  by 
Isaacs  (Ref.  [19])  is  sufficient  to  predict  lift  history  because 
the  flow  is  largely  attached  and  falls  within  the  limitations 
of  the  potential  flow  assumptions  from  which  this  theories 
was  derived.  This  is  not  the  situation  at  6  =  20°  incidence 
(Fig.  |12(b)|)  where  viscous  effects  relating  to  flow  separation 
predominate.  Consequently,  Isaacs’  model  grossly  overpre¬ 
dicts  peak  lift  for  all  cases  at  high  incidence  and  does  not  ac¬ 
count  for  the  viscous  interactions  that  lead  to  a  phase  lead  of 
the  peak  lift. 

At  6  =  20°  (Fig.  1 12(b)),  computations  predict  similar 
trends  relative  to  the  Isaacs  model  but  give  a  better  estimate 
of  the  peak  lift.  For  A  =0.8,  the  peak  lift  measured  experi¬ 
mentally  is  well-predicted  by  the  numerical  simulations,  for 
A  =  1.0  it  is  slightly  underestimated,  and  for  A  =  0.57  it  is 
largely  overestimated.  Comparing  computational  flow  fields 
and  flow  visualization  for  t/T  =  0.25  reveals  more  distinct 
large-scale  structures  as  a  result  of  the  roll  up  of  the  shear 
layer  in  the  computational  results  than  in  the  experiments.  The 
presence  of  more  coherent  vortical  structures  above  the  air¬ 
foil  leads  to  an  increased  (and  in  this  case  overpredicted)  lift. 
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Fig.  12.  Force  history  comparison  between  numerical  simulations,  experiments,  and  theory. 


Furthermore,  the  computations  do  not  predict  the  phase  lead 
of  the  experimental  lift  histories  with  regard  to  the  motion. 
With  increasing  k  or  A,  maximum  lift  no  longer  coincides  with 
maximum  dynamic  pressure  (at  t/T  =  0.25)  as  predicted  by 
the  inviscid  model,  but  occurs  earlier  in  the  period.  For  these 
cases,  the  velocity  relative  to  the  wing  is  larger  and  there¬ 
fore  the  vorticity  production  that  is  fed  into  the  shear  layer  is 
greater,  leading  to  a  steeper  rise  in  lift,  a  larger  peak  lift,  and 
stronger  viscous  interaction  in  the  shear  layer  early  in  the  ad¬ 
vancing  phase  of  the  cycle.  As  the  cycle  progresses,  stronger 
viscous  effects  might  cause  precipitation  of  deep  stall,  a  weak¬ 
ening  of  bound  circulation,  and  a  drop  in  lift  ahead  of  the  drop 
in  dynamic  pressure. 

Note  that  normalized  lift  at  0  =  20°  and  t/T  =  0  is  near  2, 
and  not  near  1.  Dynamic  pressure  at  t/T  =  0  and  t/T  =  0.5 
is  just  the  free-stream  value,  so  a  quasi-steady  lift  response  in 
this  normalization  at  those  phases  would  be  1 .  Also,  quasi¬ 
steady  lift  would  be  maximum  at  t/T  =  0.25,  where  the  in¬ 
stantaneous  dynamic  pressure  is  maximum.  In  the  computa¬ 
tion  this  is  approximately  true,  but  in  the  experiment  the  lift 
peaks  before  t/T  =  0.25,  i.e.,  at  a  lower  instantaneous  dy¬ 
namic  pressure.  It  is  conjectured  that  this  could  be  due  to 
the  partial  reattachment  of  the  suction-side  flow  shortly  after 
t/T  =  0.  As  the  blade  section  accelerates  into  the  free  stream, 
suction  is  progressively  attenuated,  and  thus  the  lift  (and  also 
the  drag;  or  more  properly,  the  normal-force)  will  peak  before 
reaching  maximum  dynamic  pressure. 

The  extent  to  which  lift  history  departs  from  the  quasi¬ 
steady  can  be  further  elucidated  by  calculating  lift  not  with 
respect  to  the  steady  free-stream  dynamic  pressure,  but  with 
respect  to  the  instantaneous  value  of  dynamic  pressure,  ob¬ 
tained  from  superimposing  the  surging  speed  of  the  blade  and 
the  free  stream  speed  of  the  tunnel.  These  results  are  pre¬ 
sented  in  Figure  [13]  In  this  formulation,  a  truly  quasi-steady 
response  would  be  identically  1  for  all  t/T,  even  in  the  case 


Fig.  13.  Lift  history  normalized  with  respect  to  the  instan¬ 
taneous  value  of  dynamic  pressure. 

ofA  =  lat^/r  =  0.75  (zero  lift  produced  at  zero  dynamic 
pressure,  or  0/0).  It  was  already  mentioned  that  at  t/T  =  0, 
the  measured  and  computed  lift  are  near  2,  and  the  measured 
lift  peaks  occur  before  t/T  =  0.25.  At  later  t/T,  lift  normal¬ 
ized  by  instantaneous  dynamic  pressure  declines  monotoni- 
cally,  especially  for  the  experimental  cases,  which  fall  below 
the  quasi-steady  value  (that  is,  below  1)  at  t/T  ~  0.6.  Ap¬ 
proaching  t/T  =  0.75,  there  is  a  rise  in  lift.  Overall,  departure 
from  quasi- steady  response  becomes  greater  with  increasing 
A.  The  lift- values  from  Isaacs’  theory  appear  artificially  large 
in  this  presentation,  because  the  Isaacs’  values  are  normalized 
by  the  static  viscous  lift  coefficient  (an  experimentally  mea¬ 
sured  force  coefficient,  not  dynamic  pressure).  Were  Isaacs’ 
formulation  normalized  by  2k a,  it  too  would  be  1  at  t/T  =  0. 
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Conclusions 

A  blade  section  in  streamwise  oscillation  was  taken  to  be  a 
model  of  a  retreating  helicopter  rotor  blade  section  includ¬ 
ing  reversed  flow  conditions.  The  lift  histories  for  differ¬ 
ent  amplitudes  of  streamwise  oscillations,  reduced  frequency, 
and  advance  ratio  were  measured  experimentally  and  com¬ 
puted  numerically  for  a  NACA0012  at  both  low  and  high  in¬ 
cidence.  Numerical  and  experimental  flow  visualization  was 
performed  to  identify  key  flow  features  and  provide  physical 
insight  into  the  lift  histories. 

For  a  blade  section  at  low  incidence,  lift  histories  for  all  of 
the  measured  cases  show  similar  trends.  During  the  first  quar¬ 
ter  of  the  semi-stroke,  lift  increases  and  reaches  a  maximum 
near  t/T  =  0.25,  when  the  dynamic  pressure  is  maximum,  in 
accordance  with  Isaacs’  theoretical  model.  At  this  incidence, 
the  flow  is  largely  attached  during  the  entire  advancing  stroke 
and  only  shows  signs  of  separation  when  the  relative  velocity 
is  lowest. 

For  a  blade  section  at  high  incidence,  lift  histories  become 
strongly  dependent  on  advance  ratio  and  reduced  frequency. 
The  peak  lift  generally  increases  with  A  and  k,  and  a  phase 
lead  arises  with  respect  to  the  motion,  likely  as  a  result  of  vis¬ 
cous  interactions.  Flow  is  fully  separated  throughout  most  of 
the  cycle.  Isaacs’  inviscid  model  strongly  overpredicts  peak 
lift  and  is  unable  to  predict  the  phase  lead.  There  is  a  fair 
agreement  between  the  computational  and  experimental  peak 
lift  for  moderate  advance  ratios  and  reduced  frequencies,  but 
computations  do  not  predict  the  phase  lead  observed  in  exper¬ 
iments. 
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